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Abstract 

Integrin signaling regulates cell migration and plays a pivotal role in developmental processes and cancer 

metastasis. Integrin signaling has been studied extensively and much data is available on pathway 
components and interactions. Yet the data is fragmented and an integrated model is missing. We use 
a rule-based modeling approach to integrate available data and test biological hypotheses regarding the 
role of talin, Dokl and PIPKI in integrin activation. The detailed biochemical characterization of integrin 
signaling provides us with measured values for most of the kinetics parameters. However, measurements 
are not fully accurate and the cellular concentrations of signaling proteins are largely unknown and 
expected to vary substantially across different cellular conditions. By sampling model behaviors over the 
physiologically realistic parameter range we find that the model exhibits only two different qualitative 
behaviours and these depend mainly on the relative protein concentrations, which offers a powerful point 
of control to the cell. Our study highlights the necessity to characterize model behavior not for a single 
parameter optimum, but to identify parameter sets that characterize different signaling modes. 
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Introduction 



Cell migration is a carefully regulated process that is essential for embryonic development and life ( Horwitz 
and Webb, 2003). As the cell moves adhesion complexes form and dissolve. Key molecules in such 



focal adhesions are integrins, large membrane-spanning molecules that bind to ligands outside the cell 



and a variety of regulatory proteins inside the cell (Arnaout et al., 2005 Banno and Ginsberg, 2008 



Campbell, 2008). Integrins are allosteric proteins that can respond to extracellular and intracellular 
stimuli and change their affinity for ligand (Hynes et al., 2002 1. The two extreme conformations, an 
open and a closed one, bind ligand with maximal and minimal affinity respectively. The extracellular 
conformational changes are accompanied by movements of the intracellular domains which lead to a 



separation of the integrin tails (Hynes et al., 2002). Binding of ligand shifts the equilibrium to the active 
open conformation. The separated integrin tails can then bind further signaling proteins and link to 
the cytoskeleton (Harburger et al., 2009). Intracellular activators such as talin and kindlins can also 



trigger integrin activation, a phenomenon that is referred to as inside-out signaling (Moser et al., 2009) 



Ligand-dependent outside-in and signaling-dependent inside-out signalling are no separate processes; 
ligand binding leads to the activation of intracellular proteins that can, in principle, feed back on integrin 
activation. In fact recent experiments show that binding of talin to the cytoplasmic tails is essential for 



ligand-dependent integrin activation (Zhang et al., 20081. In the absence of talin, interaction with ligand 



leads only to a transient activation of downstream signaling and cells fail to adhere to the substrate ( Zhang 



et al., 2008) 



Talin binds to the integrin beta-tail and stabilizes the active, open integrin conformation (Tadokoro 



et al. 



2003). Most cellular talin is unavailable for integrin binding because of self- interactions between 

These inhibitory interactions can be 
20081). PIP2 is produced by type I 



the PTB binding region and a tail region ([Goksoy et al., 2008 ) 
reheved by binding of the lipid PIP2 ([Martel et al., 2001[ [Goksoy et al 



phosphatidylinositol phosphate kinase-7661 (PIPKI) and recruitment of PIPKI to focal contacts requires 
talin binding (Paolo et al., 2002 Ling et al., 2003[ |de Pereda et al. , 2005). Ligand-bound integrins can 



stimulate the activity of PIPKI by enabling Src- mediated phosphorylation (Ling et al., 2002 Ling et al. 



2003). Src kinase binds to beta-3 integrin tails ( Arias-Salgado et al., 2005 Arias-Salgado et al., 2003) 



and ligand-dependent clustering of integrins has been suggested to trigger Src auto-transphosphorylation. 
Integrin activation may thus trigger a positive feedback loop in that activation of Src kinases and PIPKI- 
dependent talin activation and recruitment enhances integrin activation. However, the architecture of 
this positive feedback loop is further complicated by the observed competition between integrin tails 
and PIPKI for talin binding ([Barsukov et al., 2003| |de Pereda et al., 2005[ [Goksoy et al., 2008[). Thus 



Src-dependent phosphorylation of PIPKI enhances the binding of talin and PIPKI, while Src-dependent 
phosphorylation of integrin beta-tails reduces their affinity for talin and increases their affinity for other 
competing signaling protein, i.e. Dokl (Ling et al., 2003 Oxley et al., 20081. The latter effect has 



been coined integrin phosphorylation switch and has been suggested to induce a temporal switching from 
talin-dependent to Dokl-dependent integrin signaling. However, since only talin but not Dokl stabilizes 



the open, active integrin conformation (Wegener et al., 2007) it is unclear whether this switching is 
self-limiting and whether it can confer a switch in downstream signaling. The regulatory system is 
remarkably sensitive to the concentration of PIPKI: both a lower and a higher concentration impede 



talin recruitment and cell spreading (Ling et al., 2002 1. Does this help PIPKI to fulfill a dual role in first 



supporting integrin activation and then terminating integrin activation in a competition for talin ( Oxley 
et al., 2008)? Questions remain also regarding the exact mechanism as well as the purpose of these 



feedbacks. In particular, considering that ligand binding appears to be sufficient to trigger rapid and 



maximal integrin and Src activation (Zhang et al., 2008) and that PIPKI has been suggested to sequester 



talin at a later stage (Oxley et al., 2008) it is unclear why talin would be necessary for sustained integrin 
activation. 

Such questions are difficult to address by verbal reasoning alone. Mathematical modeling can help to 
integrate available isolated experimental information into a single model and permits the efficient analysis 
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and comparison of model alternatives. Given the many states and complexes that need to be considered 
the dynamics of integrin activation can be best captured by a rule-based modeling approach fFaederJ 
[eFal., 2009 ). Rule-based modeling permits the use of available information about complex protein-protein 
interactions in a precise and compact way. It is thus a convenient tool to construct a large and complex 
signaling network from a set of biochemical reaction rules (Nag et al., 2010). In decades of detailed 
biochemical analysis of proteins and their interactions most relevant rate and equilibrium constants have 
been measured. Yet, measurements are not fully accurate and the cellular concentrations of signaling 
proteins are largely unknown and expected to vary substantially across different cellular conditions. 
To account for this level of uncertainty, we devise an ensemble modeling approach to characterize the 
biologically feasible dynamic range of integrin signaling. We realize this approach by a parameter sampling 
strategy. By integrating the available biochemical information and employing an ensemble modeling 
approach we address the following questions: (i) Can our model recapitulate both modes of signaling: 
outside-in and inside-out? What are the respective parameters/control points for both signaling modes? 
(ii) What is the role and mechanism of mutual talin/PIPKI membrane recruitment during signaling? (iii) 
What is the role of integrin phosphorylation and Dokl recruitment during integrin activation? Is there 
an integrin phosphorylation switch as hypothesized previously (Oxley et al., 2008)? 



Results 

Model setup and parameterization 

Integrin activity is regulated by the interaction with the extra-cellular matrix (outside-in signaling) and 
intra-cellular cytoskeletal adaptors (inside-out signaling). We use a rule-based modeling approach, in 
order to capture the complexity of integrin signaling. This approach allows us to model integrin signaling 
compactly by a set of 30 biochemical reactions rules, which are summarized graphically in Figure [T] and 
are described in detail in the Supplementary Material. The model considers six components: integrins 
(INT), ligand (L), talin (TAL), PIP kinase (PIPKI), Src kinase (SRC) and Dokl (DOK). The rules 
describe molecular interactions, state transitions (such as protein phosphorylation), and translocation 
between cytoplasm and membrane. We consider these two compartments because recruitment of proteins 
to the membrane enhances their local concentration. By considering two compartments we can capture 
this effect without having to include space explicitly in our simulation. We do not model the production 
and turnover of proteins, as these happen mostly on a time scale different from integrin activation. Figure 
[i]A_ depicts the possible interactions (solid lines) and translocations between cytoplasmic and membrane 
compartments (dashed lines). Binding sites are specified by a circle, while interactions are depicted as 
solid lines connecting circles. Interactions which are competitive and cannot occur at the same time 
are shown as half-filled circles. Binding sites whose affinity is regulated by phosphorylation, like the 
NPxY binding motif on integrin tails, have an additional state indicated (U or P in this case). Other 
phosphorylation motifs are shown as square boxes containing either U or P while conformational states 
are also given as square boxes with either an open (O) or closed (C) state. The corresponding rule number 
which encodes the interaction or state transition is given next to each link. All state transitions which 
do not reflect binding interaction or translocations are additionally shown in Figure [TjB. Note that each 
rule encodes only the relevant biochemical context for its reaction to happen. It could therefore apply 
to many species at the same time. We have encoded the 30 reaction rules within the BNGL modeling 



language (Faeder et al., 20091 and compiled them into a system of ODEs based on the assumption of 
mass-action kinetics. This results in a reaction network of 108 species connected by 456 fluxes. 

The basic dynamics of our model are as follows. Binding of ligand and/or talin stabilizes the 
open/active conformation of integrin (rule R2, rules R3a and R4a). In case of outside-in signaling integrin 
signaling is initiated by binding of ligand (rule R2 in Figure [l] A) . Most cellular talin is unable to bind 
integrins because of inhibitory self-interactions. Talin can be activated by PIP2. Since PIP2 is highly 
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unstable, its production, decay and diffusion are not modeled explicitly. We rather assume that talin 
must be bound to membrane-bound, active PIPKI for activation (rule Rlla) and only a small fraction 
of membrane bound talin can be activated independently of PIPKI (rule Rllb). PIPKI is activated by 
SRC-mediated phosphorylation (rule R8). Src kinases are activated by trans auto-phosphorylation which 
in case of outside-in signaling is most likely mediated by juxtaposed integrin-ligand complexes. We do not 
model juxtaposition of integrins explicitly but use ligand-bound integrins as a proxy for juxta-position of 
integrins (rule R7a). The second mode of integrin activation, inside-out signaling, depends on integrin- 
independent talin activation and recruitment to the membrane. The GTPase Rapl appears to play an 
important role in the recruitment of talin (Banno and Ginsberg, 2008). Moreover, cross-talk from other 
pathways may lead to Src kinase and subsequently PIPKI activation. Since direct talin recruitment by 
Rapl/RIAM provides a straight-forward mechanism for integrin activation we focus on the more intricate 
path via Src kinase activation by cross-talk from other pathways (rule R7b in Figure [T]B). TAL, DOK 
and PIPKI can translocate between cytoplasm and membrane compartments (rule R16, R17a/b/c and 
R18) which further modulates integrin activity. 

Given the large body of experimental literature on integrin signaling measured values could be ob- 
tained for almost all rate constants (Table [T]). Nevertheless, some uncertainty remains, partly because of 
inaccurate or missing measurements, but also because of natural variations, in particular in protein con- 
centrations. Given the likely regulatory impact of different protein concentrations we sought to analyze 
the entire physiological plausible range. Thus, based on the experimental data we defined for each of the 
33 model parameters the most likely value as well as the range in which the parameter value was likely 
to lie. Here we considcrd three classes of uncertainty: parameters for which there is detailed biochemical 
data were explored over a 2-fold range, parameter values that were based on indirect or less reliable data 
were modulated over a 5-fold range and parameter values that were based on measured values for similar 
proteins were allowed to change by up to 10 fold. The maximal range to be explored was thus a 100-fold 
difference in either concentration or time-scales for one parameter. Larger changes are unlikely to be 
of relevance and may lead to effective model reductions due to scale separations. Model parameters, 
measured values and references and as well as the considered ranges which reflect parameter uncertainty 
are given in Table [l] 



Analysis of parameter uncertainty 

We use a sampling strategy to evaluate the influence of parameter uncertainty. As initial conditions we 
always used a common ground state of the network, i.e. the state that is attained without external acti- 
vators such as ligand or external Src kinase activation. We then sampled each parameter independently 
from an exponential distribution within the range specified in Table [ij This approach achieves an even 
sampling of the orders of magnitude as encoded by the fold change. In principle, our sampling space 
is very large. There are (2 x 2)^(2 x 5)"'^^(2 x 10)"'^^ ~ lO^'* possible fold-change combinations, making 
exhaustive sampling infeasible. However, we expect the space of qualitatively different model dynamics 
to be much smaller. In order to estimate the number of samples needed to obtain a comprehensive pic- 
ture of the possible dynamic range of the model for the given parameter uncertainty we monitored the 
convergence of the mean of each observable and time point in dependence of the sample size. We obtain 
convergence for 10^ parameter samples {CV < 0.1, see supplemental data). We note that this approach, 
in general, will not capture rare dynamical events. 

When we simulate ligand-dependent signaling with the 10^ different parameter combinations we no- 
tice that the level of biologically relevant output variables (i.e. open integrins, membrane-bound TAL, 
and TAL-bound integrins) varies strongly between simulations (Figure [2]). In spite of this large uncer- 
tainty we can clearly discern two time scales in the model behaviour: fast ligand-driven activation of 
integrins (occurring on the second-to-minute scale) and slow TAL membrane recruitment and TAL:INT 
and DOK:INT complex formation (occuring on the minute-to-hour time scale). However, apart from 
the time scales little can at first be said about the extent of integrin activation and its interaction with 
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talin or Dokl because 95 % of the results cover more than 50% of the possible dynamic range (Figure 
[2]B,D dotted lines). Similarly, DOK binding to integrin is highly variable when comparing its dynamic 
range (95 % of the results cover 20% of the dynamic range) to the low mean level of DOK:INT complex 
formation (less than 1%) (Figure [2] C) . 

We wondered whether the impact of parameters on the biologically relevant observables would be 
correlated. We therefore analysed pairwise scatter plots for eight observables that best capture the 
biologically interesting behaviour of the model, using the 10^ parameter samples underlying Figure [2] 
The first two observables quantify the extend of DOK and TAL binding to integrins in response to ligand 
stimulation (outside-in signaling) (FigurejsjA). The mutual dependence of TAL and PIPKI for membrane 
recruitment is captured by observables three and four (FigureJslB), while the level of integrin and PIPKI 
phosphorylation is quantified by observables 5 and 6 (Figure jsTC). The last two observables characterize 
the extend of integrin activation in response to outside-in and inside-out signaling (Figure [3]D). A formal 
definition of all eight criteria is given in the Materials and Methods section. As TAL and DOK bind to 
integrin tails in a competitive manner, we expect that TAL:INT and DOK:INT are exclusive states for 
most of the parameters sets. This is indeed confirmed by Figure[3]A, which shows beside the large marginal 
dynamic ranges of TAL:INT and DOK:INT, that both complexes are negatively correlated throughout 
the parameter samples. Positive correlation is observed between the extent to which parameter sets 
enable outside-in and inside-out signaling (Figure |3]D). Thus, whenever a parameter set achieves a high 
level of INT activation by outside-in signaling, it also allows for high levels of INT activation by SRC 
crosstalk in the absence of any ligand. 

To better understand how the eight criteria interlink, we picked the 2% extreme cases of the distri- 
bution in Figure [s] A, the parameter sets allowing either for high TALiINT formation (red points) or 
high DOK:INT formation (blue points) and marked the corresponding parameter sets also in the other 
three scatterplots. In this way we find that TAL:INT complex formation and the ability for high levels 
of outside-in and inside-out signaling is strongly linked (compare red points Figure [s] A and D). Also 
the level of INT phosphorylation is lower for parameters which achieve high levels of INT activation 
(i.e. open integrins) (Figure |3]C), while the level of PIPKI activation is high (Figure [3]B). Interestingly, 
parameters which favor DOK:INT complex formation do not allow for inside-out signaling (compare blue 
points Figure [3] A and D). INT phosphorylation in this case is also limited but significantly higher than 
for the TAL:INT favoring parameter sets (FigurejsjD). The distribution of PIPKI activity is rather broad 
(Figure [3] B). Also the membrane recruitment of PIPKI is in most cases independent of TAL as seen by 
the distribution of blue points in the lower part of the vertical axis in Figure |3]B. Parameters which allow 
for high levels of INT activation have a clear tendency for TAL-dependent PIPKI recruitment in this 
case. Thus, it seems that for all eight criteria we can isolate two sets of strongly correlated parameters 
which can be linked to two distinct dynamical regimes of the model: 

Group 1 is characterized by a strong potential to support INT activation either by outside-in or inside- 
out signaling, high levels of TAL:INT complex formation, a strong dependence of PIPKI recruitment 
on TAL, low levels of INT phosphorylation and high levels of PIPKI activation. 

Group 2 is characterized by the inability to support inside-out signaling, high levels of DOK:INT com- 
plex formation, no dependency of PIPKI recruitment on the presence of TAL and higher levels of 
INT phosphorylation compared with group 1. 

We next sought to define the parameter ranges that would correspond to the two distinct model 
behaviours and identify those parameters that would affect the two model behaviours the most. To this 
end we selected the parameter samples that corresponded to each group and compared the distribution 
of the single parameters (Figure |3]E, group 1 (red boxes) or group 2 (blue boxes)) with their original 
sampling distribution (Figure [sjE, grey boxes). We reasoned that sensitive parameters would be sampled 
from a strongly restricted range. The extent of the deviation was computed as the maximal difference 
between the cumulative distribution function (CDF) of a parameter and the CDF of its uniform sampling 
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distribution. The absolute maximal difference, also called the Kolmogorov-Smirnov test statistics, was 
used to rank parameters for their relative influence on the group behavior (Figure [3]F, red bars for group 
1, blue bars for group 2, both versus the uniform sampling distribution, and gray bars for the difference 
between group 1 and group 2). The strongest deviation from the uniform sampling distribution is seen 
in the parameters for the relative total protein concentrations. Group 1, which is characterized by the 
capability of inside-out signaling, has a clear preference for higher TAL levels compared with the levels 
of integrin, DOK and PIPKI. Additionally, the PIPKI-dependent TAL activation rate is enhanced, while 
the PIPKTindependent activation is diminished and the binding and unbinding of TAL and PIPKI is 
shifted to faster time scales. Together this indicates a crucial role for TAL and PIPKI for the behavior 
of group 1. The same parameters tend to be distributed differently in group 2. Here the distribution of 
the total concentration ratios is just the opposite with a lower TAL concentration in comparison to INT, 
DOK and PIPKI levels. Some of the kinetic parameters are of characteristic importance, too. Thus the 
time scales for the binding and unbinding of DOK to phosphorylated INT are much faster on average, 
further disfavouring the formation of TALiINT complexes and allowing for higher levels of DOK to INT 
binding. Interestingly the time scales for ligand-integrin interaction are lower on average, suggesting that 
a more stable ligand INT interaction correlates with higher extents of DOK binding. This result can be 
understood in light of the particular biochemical constraints: ligands need to stabilize open INT for DOK 



to bind to integrin tails because DOK cannot stabilize INT in its open conformation (Wegener et al. 



2007). Overall this analysis highlights that the qualitative behaviour of biological networks depends 
only on few parameters. Interestingly, the two different dynamic regimes uncovered by our analysis 
(group 1 and 2 parameter sets) are mainly deflned by the relative total concentrations of TAL, DOK 
and PIPKI. Cells have powerful mechanisms in place to adjust relative protein concentrations. These 
thus provide excellent control points for cells to define its cellular dynamics. We thus hypothesize that 
integrin signaling operates in both regimes by altering the relative levels of TAL, DOK and PIPKI. In 
the following we investigate their biological role in more detail. 

Integrin phosphorylation switch 

Detailed biochemical measurements have led to the proposition of an integrin phosphorylation switch 



(Oxley et al., 2008). Thus measurements revealed that Src-dependent phosphorylation of integrin beta- 



tails lowers the integrin-talin affinity by about 2-10-fold (Ling et al., 2003 Oxley et al., 2008) while 



enhancing the affinity of Dok-1 for the same binding site some 400- fold (Oxley et al., 2008). Some 



experiments also reported an increased affinity of PIPKI for talin once PIPKI had been phosphorylated 



by Src kinases (Ling et al., 2003 de Pereda et al., 2005). Src-dependent phosphorylation of integrin 



tails and PIPKI upon integrin activation was thus suggested to result in a time-dependent exchange of 
TAL for DOK on integrin tails and the sequestering of TAL in complexes with PIPKI, coined integrin 



phosphorylation switch (Oxley et al., 2008). We use our simulation to explore this hypothesis for the two 
groups of parameter sets identified above. Thus Figure[4]A-D shows the extent of integrin phosphorylation 
and binding to DOK or TAL over time and dependent on the DOK/TAL ratio for one set of parameters 
characteristic for each group of parameter sets. The characteristic parameter set represents the mean 



value of all parameter samples belonging to each group Table [ST] As expected, the level of TAL:INT 
exceeds the level of DOKTNT in group 1 (Figure [3] A,B), and vice versa for group 2 ( Figure[3]C,D). The 
level of integrin phosphorylation is higher for group 2 because in case of group 1 parameters TAL binding 
shields the NPxY phosphorylation motif in integrin tails and thus prevents phosphorylation, while in 
case of group 2 parameters the unphosporylated NPxY motif is accessible and bound DOK protects the 
phosphorylated NPxY motif from dephosphorylation. In spite of strong integrin phosporylation for group 
2 parameters we, however, do not observe a temporal switching from TAL:INT to DOK:INT association 
(Figure |3]C). Instead, both complexes coexist for comparable levels of TAL and DOK in case of group 2 
and balance in favor of TALTNT in case of group 1 (Figure [s] A, C). 

Although there is no temporal phosphorylation switch, its is possible to switch between TAL and 
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DOK bound forms of integrin by regulating the ratio of total concentrations of DOK and TAL (Figure |4] 
B and D). For group 2 the switching point is reached at equal DOK and TAL total concentrations which 
is close to the likely physiological protein concentrations (dotted vertical lines). For group 1 parameters 
on the other hand the switching point is reached at 1000-fold higher levels of DOK over TAL which is 
physiologically unrealistic. Thus, in group 1, TALiINT complex formation seems robust to the DOK-to- 
TAL ratio, while for the parameter sets in group 2 physiological changes in the protein concentrations 
could lead to a switch. Figure|4]B and D confirm our previous observation that the level of phosphorylated 
integrin negatively correlates with the total level of TAL and DOK-bound integrins. 

While we can exclude a temporal switch in DOK and TAL binding of integrin tails we note a temporal 
switch in integrin phosphorylation for group 1 parameters (Figure |4] A) . Thus as talin accumulates at 
the membrane and binds to integrin tails the extent of integrin phosphorylation decreases after having 
reached an early maximum of about 60% of phosphorylation. Integrin signaling is thus transient. To 
see whether this is a general behavior for group 1 parameters and wether its also hold for group 2 
parameters we analysed the entire parameter sets belonging to either group ( Figure |4]E). We used two 
measures for transient responses: sensitivity towards stimulation and the precision with which the extent 
of integrin activation (E), INT phosphorylation (F) and TAL:INT complex formation (G) return to pre- 
stimulation levels (see materials and methods section for details). Transient responses are characterized 
by high sensitivity and precision while sustained responses are characterized by high sensitivity and low 
precision. A low sensitivity would characterize unresponsive systems. We observe a clear separation for 
the two parameter sets in the scatter plots in Figure [4] E-G. Thus group 1 parameter sets result in self- 
limited integrin phosphorylation; phosphorylation levels, however, do not return to base line but assume 
some intermediate level (i.e. 50% phosphorylation. Figure |4] A-l-F). Group 2 parameters on the other 
hand enable a very sensitive, sustained phosphorylation response. Even though integrin phosphorylation 
is necessary for downstream signaling, integrin activation is typically monitored as the extent to which 
integrins assume the open conformation. Interestingly, this alternative measure of integrin activation 
gives very different results (Figure |4]E). Thus group 1 parameters lead to strong and sustained integrin 
activation (even though phosphorylation is self-limiting) while Group 2 parameters lead to intermediate (if 
sustained) integrin activation. In fact all simulated 10^ parameter combinations result in sustained levels 
of open, active integrins as indicated by the distribution of all parameter sets along the diagonal (Figure 
|4]E). As expected from the previous analysis, the establishment and temporal dynamics of TAL:INT 
signaling correlates well with integrin activity, showing high sensitivity and sustained response for group 1 
parameters (Figure|4]G). Phosphorylation and downstream signaling can thus be independently controlled 
from the conformational change that initiates integrin activation. 



Mutual regulation of TAL and PIPKI membrane recruitment 

TAL and PIPKI are cytoplasmic proteins that need to be enriched at the membrane to serve their 
function in integrin activation. However, it has been difficult to explain how they become enriched at 
the membrane as both bind to membrane lipids at most weakly. In the following we focus our analysis 
on group 1 parameters since only these permit strong INT:TAL complex formation (Figure [s] A). Figure 
[3]B suggested that PIPKI and talin mutually enhance each others membrane recruitment, and indeed 
for very low levels of either TAL or PIPKI there is no enhanced membrane recruitment (Figure [s] A and 
B). TAL enhances PIPKI recruitment by acting as a membrane adaptor that binds both the membrane 



and PIPKI. PIPKI, on the other hand, enhances TAL recruitment in the model because as a dimer ( |Rao 



et al., 1998) it multimerizes TAL and thus enhances its affinity for the membrane. Moreover, it activates 
TAL by producing PIP2 and thus enables its binding to integrin tails. Since the PIPKLTAL interaction 
does not depend on integrin activation both accumulate at the membrane to a considerable extent also in 
the absence of ligand (Figure [5]G). Integrin activation then enhances recruitment, mainly by permitting 
INT:TAL complex formation. As observed previously these complexes only form minutes to hours after 



ligand addition even though open integrins emerge within seconds (Ling et al., 2002) 
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The extent of integrin activation strongly depends on the total talin and PIPKI concentrations (Figure 
[5] A and B). Thus the fraction of active integrins increases as the talin concentration is increased and 
reaches near maximal levels at the likely physiological concentration of 10 /iM (Figure [5] A). Both the 
extent of integrin activation and talin recruitment to the membrane exhibit a bell-shaped dependency 
on the PIPKI concentration (Figure [s] B). At low PIPKI concentrations too little scaffold is available 
for talin to be multimerized at the membrane while at very high concentrations all talin becomes se- 
questered in complexes with PIPKI and is no longer available for integrin activation. Similar as discussed 
for TAL/DOK there is thus no temporal but a concentration-dependent switching point for integrin ac- 
tivation. The bell-shaped model prediction is in agreement with the experimental observations: mild 
over-expression of PIPKI enhances talin recruitment while strong over-expression of PIPKI inhibits talin 
recruitment (Ling et al., 20021 and leads to the dissolution of focal adhesions (Paolo et al., 2002 1. This 
experimental observation raised the question as to why physiological levels are suboptimal for talin re- 
cruitment. Based on our simulations we note that the curves for talin recruitment and integrin activation 
do not coincide. Maximal integrin activation is achieved for lower levels of PIPKI compared with the 
levels for maximal TAL recruitment. The mean value of group 1 parameters (which were selected for 
highest TAL:INT formation) coincides with maximal integrin activation (1 PIPKI, Figure [5]B, dot- 
ted line). We thus propose that the signaling network was optimized to maximize integrin activation in 
response to ligand binding, and that this maximum is achieved when the extent of talin recruitment is 
sub-optimal as observed in experiments. 



Outside-in and inside-out signaling 

Integrins are unusual receptors in that they can be activated both by external signal (outside-in) and 
internal signals (inside-out signaling) . As discussed above both modes of signaling are observed only with 
group 1 parameters (Figure |3]D) and we will therefore now focus on these parameter sets. Both signaling 
modes result in a similar extent of integrin activation yet the kinetics are very different. Thus outside-in 
signaling follows a biphasic time course with an early talin-independent phase and a later talin-dependent 
phase (Figure|6]A, compare black solid and dashed lines) while inside-out signaling exhibits only the later 
talin-dependent signaling phase (Figure [6] A, compare magenta solid and dashed lines). The increased 
talin dependency of inside-out signaling becomes apparent also in Figure [6] B where we record steady 
state integrin activation in dependence of total talin levels for both signaling modes. 

Figure |6] C illustrates the extent to which talin supports outside-in signaling at low ligand densities. 
High levels of open integrins appear already for low ligand levels but are absent in the TAL knock-out. 
As expected, SRC levels play a similar role for inside-out signaling as ligand levels for outside-in signaling 
(Figure [6] D) since total SRC levels are a proxy for the amount of SRC which becomes activated by cross- 
talk signaling. Both INT activation curves are qualitatively similar, except that full INT activation is only 
achieved in the presence of high ligand concentrations. Interestingly, both curves show an intermediate 
range of high integrin activation {L/Kjj ^ 10~^/iM, Figure[6]C) and for inside-out signaling lower levels 
of SRC achieve highest INT activation (SRC - lO^^M, Figure |6]D). In both cases, membrane levels 
of PIPKI increase with the amount of INT activation. Since SRC binds and sequesters PIPKI at the 
membrane, it indirectly attenuates the role of PIPKI in recruiting TAL. Thus, as observed in Figure [5] 
B, there is an optimal level of PIPKI at the membrane for INT activation, which is lower and depends 
on the level of ligand and SRC. Despite this intermediate maximum, the model predicts high levels of 
integrin activation over a wide range of ligand and Src kinase concentrations. Thus, while PIPKI levels 
appear to be optimized for integrin activation, the signaling mechanism appears to be much more robust 
to variations in ligand and SRC levels, so that integrin signaling happens almost in a none-or-all fashion 
(Figure [6] C, D). 
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Discussion 

We have developed a rule-based model to explore the roles of talin, PIPKI, and Dokl in integrin activation 
and signaling. In spite of much detailed biochemical data a significant level of uncertainty remains 
because measurements are not fully accurate and the concentrations of signaling proteins are largely 
unknown and expected to vary substantially across different cellular conditions. It is generally infeasible 
to sample the entire physiological parameter space. We therefore addressed this challenge in modelling 
signaling pathways by restricting the sampled parameter space according to the quality of the available 
measurements. The width of the ranges was based on general rules as to the accuracy of different 
experimental methods. This restriction enabled us to systematically analyse model dynamics with only 
10^ samples of parameter sets. 

Interestingly, in spite of the wide sampling we observed only two different qualitative behaviours: 
one group of parameter sets supported integrin activation both by outside-in or inside-out signaling, 
and was characterized by high levels of talin-integrin complex formation, a strong dependence of PIPKI 
recruitment on talin, low levels of integrin phosphorylation and high levels of PIPKI activation. The 
other group of parameter sets did not support inside-out signaling, and was characterized by high levels 
of Dokl-integrin complex formation, no dependency of PIPKI recruitment on the presence of talin and 
higher levels of integrin phosphorylation. 

Importantly, the most decisive parameters were the relative cellular protein concentrations. Accord- 
ingly, the total concentrations of network components can serve as powerful control points to achieve 
distinct network dynamics. Thus for low talin concentrations we predict impaired inside-out signaling, 
yet much Dokl-integrin complexes form in response to outside-in signaling which will direct downstream 
signaling via recruitment of a range of proteins, including Ras-GAP, the adaptor protein Nek, the non- 



receptor tyrosine kinase Csk, and the phosphatase SHP2 (Ling et al., 2005). In the presence of high 
talin concentrations cells can trigger both outside-in and inside-out integrin signaling and talin can link 
integrins with the cytoskeleton. Our modeling suggests that under physiological outside-in signaling con- 
ditions both signaling platforms co-exist and that the equilibrium is mainly determined by the talin/Dok 
and talin/PIPKI ratios. 

Previous experimental observations led to the suggestion of an integrin phosphorylation switch. Bio- 
chemical measurements revealed that phosphorylation of integrins and PIPKI would lower the talin- 



integrin affinity and enhance PIPKI-talin binding (Barsukov et al., 2003 de Pereda et al., 2005 Goksoy 



et al., 2008|) . Moreover, the Dokl-integrin affinity was enhanced some 400-fold (Ling et al., 2003 Oxley 



et al., 2008). Based on our parameter screen we can dismiss the idea of a temporal phosphorylation 
switch. However, we note that a switch can be achieved if relative protein concentrations are changed. 
Such change in protein concentration must be the result of protein expression rather than membrane re- 
cruitment as the signaling induced recruitment dynamics were included in the model and are insufficient 
to trigger the switch. Membrane recruitment of talin and PIPKI itself poses an interesting conundrum in 
that both appear to depend on each other, yet both bind membrane lipids at most weakly. Based on our 
simulations we propose that PIPKI enhances talin recruitment to the membrane by offering a dimeric 



scaffold (Rao et al., 1998) that multimerizes talin and thus enhances the effective talin- membrane affinity. 
Talin, on the other hand, recruits PIPKI to the membrane by binding to the scaffold. The concentration 
of PIPKI needs to be finely balanced for successful integrin activation. Thus, at low PIPKI concentrations 
there is insufficient scaffold for membrane recruitment and insufficient PIP2 production for efficient talin 
activation while at very high PIPKI concentrations active talin will be sequestered by the kinase and 
will not be available for binding to integrin tails. Therefore, the PIPKI level is one of the main control 
points for talin and integrin activation. Slight overexpression of PIPKI enhances talin recruitment in 



experiments (Ling et al., 2002). While this suggested that the interactions are not optimized for integrin 
activation, the simulation now reveals that maximal integrin activation is achieved already at a lower 
talin concentration, and the protein network thus appears to be optimized for integrin activation rather 
than talin recruitment. 
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The simulation further revealed that the opening of the integrin conformation not necessarily leads 
to integrin phosphorylation. Because binding shields the integrin phosphorylation motif, strong recruit- 
ment of either Dokl or talin to integrin tails will counteract integrin phosphorylation. Therefore integrin 
phosphorylation is indirectly regulated by recruitment of talin and Dokl. We observe that conditions 
which enable high levels of Dokl signaling also allow for a sensitive, sustained phosphorylation response 
while talin-dependent signaling allows only for transient levels of high phosphorylation. Since talin is the 
more stable signaling platform, its stable interaction with integrin tails does not allow for high levels of 
integrin phosphorylation but a sustained level of open integrins. Our modeling therefore suggest, that 
integrin activation and integrin phosphorylation serve as independent nodes in integrin signal propaga- 
tion. Finally our simulation suggests that outside-in and inside-out signaling differ in the kinetics of the 
activation process. Thus while outside-in signaling is biphasic with a rapid talin-indcpendent and a slow 
talin-dependent activation phase, there is only the late talin-dependent activation phase in inside-out 
signaling. 

It has long been recognized that model behaviours should be analysed not for a single measured 
parameter set alone, but over a wider parameter range. Ensemble modeling strategies combined with 
various clustering algorithms have been used to reveal parameter sets that yield similar qualitative results 

The network studied here is unusual in that almost all kinetic 



(Melke et al., 2006 Kuepfer et al., 2007 



parameters have been measured (Table 



In. However, in agreement with previous studies we notice that 



the qualitative model behaviour is not particular sensitive to most kinetic parameters, somethings that 
has been referred to as sloppiness (Gutenkunst et al., 20071. More important than the kinetic parameters 
are the protein concentrations which can differ between cells and conditions and which may give rise to 
heterogeneous functionality on the single cell level (Spencer et al., 2009). The large impact of protein 
concentrations provides cells with a powerful level of control and enables cells to assume to different 
modes of signaling in the context of our model here. 

Overall, our study illustrates how detailed biochemical measurements together with a thorough com- 
putational analysis can be used to reveal the qualitative behaviours of biological signaling networks. 
Without biological data the space of possible parameter combinations is too large to be explored. Yet 
without the computational model it is impossible to integrate biological knowledge to grasp the dynamic 
range of biological networks and to estimate the possible impact of parameter variations under different 
cellular conditions. 



Materials and Methods 

Model setup We use a rule-based approach to explore the roles of talin, Dokl, and PIPKI in integrin 



signaling dynamics (Faeder et al., 2009). Rule based modeling associates a pattern with a given bio- 
chemical reaction. The pattern specifies all relevant species for which the reaction is applicable. By that, 
rule based modeling gives a compact description of a biochemical reaction system and avoids redundant 
information, which may arise by the combinatorial complexity of most signaling systems. Our integrin 
signaling model consists of 6 molecule types and 29 reactions rules which were translated into a set of 108 



species and 454 individual reactions using BioNetGenerator (Faeder et al., 2009). There are two reaction 
compartments: plasma membrane and cytoplasm. The membership to a certain reaction compartment 
was specified by molecule states and the exchange fluxes between compartments were adjusted on a per 
species basis in order to account for different compartment volumes (see Supplement for details). See 
Supplemental Data for a complete model description and the set of reaction rules. Time course and 
steady state solutions of the set of ODEs were computed numerically with the odel5s integrator of MAT- 
LAB. The ground state (GS) is defined as the steady state solution in the absence of any ligand or SRC 
cross-talk activation. It served as the initial condition for all simulations of integrin activation. Outside-in 
(01) signaling was simulated by setting Ltot = 15, while inside-out (10) signaling was achieved by setting 
Ltot — and switching on rule 7b which leads to an integrin independent SRC activation. 
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Parameterization and treatment of uncertainty We collected parameter values for all biochemical 
reactions from the literature and parameterized the model according to Table [l] (see Supplement for 
details). Due to the heterogeneity of the biochemical data/essays and intrinsic uncertainty in reaction 
rate constants, we additionally defined a biologically plausible range for each rate constant. We sampled 
parameters within the ranges defined by the fold changes given in Table [l] in order to characterize the 
possible model behavior within the range of parameter uncertainty. Parameters were sampled uniformly 
on a logj^Q scale within the range kg/ f < k < kof, where / is the fold change and fep the basal parameter 
value, both stated in Table [T] In order to evaluate the number of parameter samples needed to cover 
the dynamical range of the model, we checked the convergence of the time-dependent mean of all model 
observables by a blocking procedure. Convergence was measured by the coefficient of variation and was 
achieved for approximately 10^ samples {CV < 0.1). 



Definition of output criteria We use the following eight criteria to characterize model behavior. 

1. DOK:INT formation in response to stimulus: 

Ci = {DOK:INToi-DOK:INTGs)/INTtot 

2. TAL:INT formation in response to stimulus: 

C2 = {TAL:INToi-TAL:INTGs)/INTtot 

3. PIPKI-dependent TAL recruitment: 

C3 = {TAL'Sr -TAL'j^:';;:,)lTALt,t 

4. TAL-dependent PIPKI recruitment: 

d = {PIPKI^'f"" ~ PIPKI^^'i"^)/PIPKItot 

5. PIPKI phosphorylation in response to stimulus: 

Cg = {PIPKI^\'°' ~ PIPKI^f')/PIPKItot 

6. INT phosphorylation in response to stimulus: 

Ce = {INTS^°' -INT^^''')/INTtot 

7. Integrin activation by outside-in signaling: 

Cr = {INToi-INTGs)IINTtot 

8. Integrin activation by inside-out signaling: 

Cs = {INTio-INTGs)/INTtot 

Here, Pipki and Tal refers to the steady state of the PIPKI or TAL knock-out after ligand application. 
Mem indicates the membrane bound forms and Phos the phosphorylated forms. We additionally defined 
the sensitivity S and the precision P of an observable Y as 

S=\{Y^'-- ^YGs)/Yt,t\ 

and 

P=\{Yoi-YGs)/Ytotr\ 
Here, Max refers to the maximum of Y{t) for < i < 4h. 
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Characterization of parameter classes Based on the first two criteria, we selected the 2000 pa- 
rameter sets (2%) which achieved the highest levels of TAL:INT (group 1) and DOK:INT (group 2). 
In order to identify the relevant parameters which characterize both groups, we further compared the 
distribution of group 1 and group 2 parameters set with their corresponding uniform sampling distribu- 
tions and among each other. We computed the empirical cumulative distribution function {CDF{x)) of 
parameter p for the sampling distribution and in each of the two group distributions. We rank parameters 
for their influence on group behavior based on the maximal deviation between all three CDFs (i.e., the 
Kolmogorov-Smirnov test statistics D). 

CDFP{x)\ 
CDFP{x)\ 
CDFl{x)\ 

Here subscripts specify the group (51,^2) and sampling distributions (s). 
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Figure Legends 



Figure 1. A rule-based model of integrin signaling. (A) Contact and localization map. (B) 
Additional state transition rules. The model considers six molecules: ligand (L), integrin (INT), Dok 
(DOK), talin (TAL), Src kinase (SRC) and pip kinase (PIPKI). Each molecule has a set of binding site 
and a set of possible states. Reactions take place within three compartments denoted as extracellular 
domain (EC), plasma membrane (PM) and cytoplasm (CP). DOK, TAL and PIPKI can translocate 
between the CP and PM compartments, while INT and SRC only reside within the PM compartment. 
L is extracellular. Possible binding reactions (annotated by the rule number) and involved binding site 
are indicated. Note, that some binding sites are competitive, e.g., the NPxY motif on INT can be either 
DOK or TAL bound. For a formal definition of the reaction rules see Supplemental Data. 



Figure 2. Integrin signaling dynamics within biologically feasible parameter ranges. 

Average time courses of integrin activation (A), TAL membrane recruitment (B), DOK bound INT (C) 
and TAL bound INT (D). Median (black fine), 50 % (dashed line) and 95 % (dotted line) data intervals. 
To better visualize the temporal order of integrin activation, time is represented on a log-scale. 
Parameters are varied according to the ranges given in Table [T] All simulations were performed with 
L/Kd = 50. Results are given for 10^ samples. 



Figure 3. Two groups of qualitatively different model behavior. (A-D) Scatter plots of eight 
criteria used to evaluate model behavior. (A) Ligand-induced DOK:INT versus TAL:INT complex 
formation. The 2% extreme ends of the distribution are colored either red (TAL:INT) or blue 
(DOK:INT) and reflect two qualitatively different model behaviors. Both groups of samples are mapped 
into panel B-D to visualize the dependencies between criteria. (E) Box plots of the parameter samples 
of each group. Ranges of the uniform sampling distributions, as stated in Table [T] are indicated by grey 
boxes. (F) Maximal difference in the cumulative distributions between sampling parameters and 
selected parameter sets: TAL:INT (blue), DOK:INT (red) and the difference between both groups 
(grey). 



Figure 4. Integrin phosphorylation and complex formation for average group 1 and group 
2 parameter values. Both groups are colored according to Figure |3] red for group 1, (panel A and B) 
and blue for group 2 (panel A and C). (A,C) Time course of integrin phosphorylation as well as 
DOK:INT, TAL:INT and TAL:PIPKI complex formation. (B,D) Steady state levels of the 
corresponding species as a function of the ratio of DOK/TAL. Vertical line: mean DOK/TAL ratio in 
both groups. (E-G) Sensitivity and precision, as defined in Material and Methods section, of open INT 
(E), phosphorylated INT (F) and TAL:INT complex (G). 
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Figure 5. Membrane recruitment of TAL and PIPKI through ligand induced integrin 
activation. (A) Fraction of membrane-recruited PIPKI and open INT as a function of total TAL levels. 
(B) Fraction of membrane-recruited TAL and open INT as a function of total PIPKI levels. (C) Time 

course of membrane recruitment of TAL. PIPKI and TAL/PIPKI complex. All results are based on the 
parameter set characterizing group 1 behavior. Vertical lines: group 1 TAL (A) and PIPKI (B) levels. 



Figure 6. Integrin activation by outside-in and inside-out signaling. (A) Time course of INT 
activation. Outside-in (01) signaling (black), inside-out (10) signaling (magenta), wild type (sohd), 
TAL knock out (dashed). (B) Fraction of open INT in dependence of total TAL. (C) Fraction of open 
INT in dependence of L/Kd- (D) Integrin activation by inside-out signaling as a function of SRC 
levels. Vertical line: mean group 1 parameter values. 
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Supplementary information 



Rule-based model: molecule types and total concentrations 

Ligand The concentration of the hgand molecule defines the density of integrin binding sites. However, 
the decisive quantity in our model is the binding site density relative to the binding site affinity: Ltot / Kd . 
We have explored ranges between 10~^'^ < Ltot/Ku < 10'' and used a value of 50 for all time course 
simulations. Given a Kd = 0.3/iM, this would correspond to ligand concentration of Ltot = 15/iM, which 



is among experimentally used ligand densities (Vitte et al., 2004 Cox et aZ.(2001) 



L(I) 



single integrin binding site 



Integ rin There are typically about 5 x 10^ integ rins on a fibroblast (Neff et al., 1982 Akiyama et al., 
1985) and about 1 — 30x 10^ integrins o n CHO cells (Wiseman et al., 2004 1. The surface area of CHO cells 
has been measured as 1000-2500 /im^ (Wiseman et al., 2004 1. Based on these measurements we expect 
a density of about 40-3000 integrins per /xm^ and we shall use 500 integrins per ^ir? as this corresponds 
best to both the measurements in fibroblast and CHO cells. If we assume 20 nm as the integrin-ligand 
binding distance then a density of 500 integrins per /im^ corresponds to a concentration of about 40 //M. 



INT(Conf-0~C,L,Sl,NPXY-U~P) 



Conf- 

L 

SI 

NPXY- 



0-C 



U~P 



Open and closed integrin conformation 
ligand binding domain 

SRC/FAK beta3 cytoplasmic tail binding domain 

Talin/Dok binding domain (Calderwood et al., 20021; gets phosphorylated by 
active Src (Anthis et aL(2009) 



Dokl 



DOKl(PTB) 

PTB Integrin binding domain 

Talinl 



TAL1(F2,F3- 
F2 

F3-0~C 



0~C) 

binding domain for PIP2 (membrane recruitment of talin) (Saltel et a/. (2009). 
PTB domain, bind integrin tails and PIPKI, is shielded in closed conformation 



(autoinhibition) (Ling et al., 2003 



de Pereda et al, 2005 



Moser et al., 2009) 



Src kinase We have shown previously that the experimentally observed kinetics of Src kinase activa- 
tion and deactivation can best be reconciled with a Src density of 300 per iivn? which corresponds to a 
Src kinase concentration of about 25 /iM and if Src kinases can interact with all integrin conformations 
(Felizzi and Iber, submitted). 



SRC (SH2 , SH3 , Y419~U-P) 
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SH2 interaction domain with CSK 

SH3~U~P domain can directly interact with betaS cytoplasmic tails 
PIPKI7 

PIPKI(SPLH,Y644~U~P) 



Y644-U-P phosphorylted by SRC (Ling et al, 2003) 

SPLH binding domain to F3 domain of talin, Y644 regulated 

Rule-based model: reaction rules, parameter values and fold changes 

Rule la: Integrin opening and closure Integrins are allosteric proteins that exist in many confor- 



mations ( Hynes et al., 2002 j . For simplicity we only consider the two extreme conformations, a closed. 



inactive conformation and an open, active conformation. The conformational equilibrium constant Kc 
can be derived from the observation that in the absence of ligand only about 2 — 3% of all integrins 



are in the open, high affinity conformation ( Tadokoro et al., 2003j ). The underlying transition rates 



have not been measured but a reasonable time scale is 1/10 — 1 sec. INT(Conf ~C,L,NPXY~?) <-> 
INT(Conf~0,L,NPXY-?) klon, kloff 

Rule lb: Integrin closure Applies to all open-DOK-bound-integrins irrespective of the phosphory- 
lation state. 

INT(Conf-0,L,NPXY-?!l) .D0K(BD!1) -> INT(Conf ~C,L,NPXY~?) + DOK(BD) kloff 

Rule 2: Ligand binding The dissociation constants for ayPa integrins in the open conformation has 
been measured as Kd ~ 0.3/LtM (Plow & Ginsberg(1981; Faull & Ginsberg(1995). As discussed in detail 



elsewhere (Iber and Campbell, 2006) the affinity of integrins for surface-bound ligand is much lower than 
for hgand in solution (Kuo & Lauffenburger(1993; Moy et al. {1999; Bell(1978). Accordingly we use 500- 
fold lower affinity constant. For such a lower affinity ligand binding is limited, even at the high ligand 
densities that are employed in experiments. We used fco„ — l^M^^ s~^ as on-rate and we set the off-rate 
according to the dissociation constant. 

INT(Conf-0,L) + L(I) <-> INT(Conf ~0,L! 1) .L(I ! 1) k2on, k2off 

Rule 3a: Interaction of unphosphorylated integrin and talin Talin binds to unphosphorylated 
integrin (3 tails with at a very low rate, kZaon = 5 x 10'^ s~^. The off-rate was measured as 



/c3ao// = 5 X lO"-^ s-i so that Kd = O.l/^M ( jCalderwood et al., 2002[ ). 
TAL(Loc-PM,Conf ~D,BD) + INT(Conf ~Q ,NPXY-U) <-> 
TAL(Loc~PM,Conf~D,BD!l).INT(Conf~0,NPXY-U!l) kSaon, kSaoff 

Rule 3b: Interaction of unphosphorylated integrin and Dok 
DOK(Loc~PM,BD) + INT(Conf -D , NPXY~U) <-> 
DOK(Loc~PM,BD!l) .INT(Conf~0,NPXY~U!l) k3bon, kSboff 

Rule 4a: Interaction of phosphorylated integrin and talin Phosphorylation of the integrin beta 



tail reduces the affinity of binding 2- fold (Oxley et al., 2008). 
TAL(Loc-PM,Conf-0,BD) + INT(Conf ~0 ,NPXY-P) <-> 
TAL(Loc-PM,Conf~D,BD!l).INT(Conf~0,NPXY-P!l) kSaon, 2*k3aoff 
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Rule 4b: Interaction of phosphorylated integrin and Dok (Oxley et al., 2008) 
DOK(Loc~PM,BD) + INT (Conf ~D , NPXY~P) <-> 
DOK(Loc-PM,BD!l) .INT(Conf~D,NPXY~P!l) kSbon, k3boff/400 



Rule 5: Src opening and closure In the resting state Src kinases are predominantly inactive (90 — 
95%) (Roskoski, 2004). To obtain such a conformational equilibrium we set our rates for the opening and 
clamping to fcSon = ls~^ and fc5o// = lOOs^^ respectively. 
SRC(Conf-C,SH3,Y419~U) <-> SRC (Conf -0,SH3,Y419~U) kSon, kBoff 



Rule 6a: Src integrin interaction The concentration at which half-maximal binding of Src kinases 
and integrins is achieved has been determined as EC^q ^ 5 — lO/iM ( Arias-Salgado et al., 2003) and 
accordingly we use as dissociation constant — 7.5/iM for the Src kinase-integrin interaction. 
SRC ( Conf -0,SH3) + INT(BD) <-> SRC(Conf ~0,SH3! 1) . INT(BD! 1) k6on, kSoff 



Rule 6b: Src PIPKI interaction 

SRC ( Conf ~0,SH3) + PIPKI(Loc~PM,Sl) <-> SRC(Conf~0 ,SH3 ! 1) . PIPKI (Loc~PM, SI ! 1) k6on, kSoff 



Rule 7a: integrin-dependent Src phosphorylation Activation of Src kinases is achieved by trans 
autophosphorylation on Tyr-418 (Harrison(2003) and thus requires juxtaposition of open Src kinases, 
presumably by binding to ligand-bound integrins ( [Arias-Salgado et al., 200*3 ). We do not model juxta- 
position of integrins explicitly, but require ligand bound integrins. 
SRC(SH3!1,Y419~U) .INT(BD!1,L!+) -> SRC(SH3 ! 1 , Y419~P) . INT(BD ! 1 ,L ! +) k7 



Rule 7b: cross-talk dependent Src phosphorylation Inside-out signaling is simulated by setting 
lOFLAG^l and Ltot = 0. 

SRC(SH3! 1,SH3,Y419-U) -> SRC(SH3! 1,SH3,Y419~P) I0FLAG*k7 



Rule 8: PIPKI phosphorylation The K„i for the opti mal Src kinase substrate has been determined 
as K,n = 30^M, and accordingly we use kcat = 10 s"-"^ (Songyang et al., 1995). In the absence of 
experiments that would suggest otherwise, we choose the same Km and kcat values for all Src-mediated 
reactions, i.e. also for the phosphorylation of integrin /3 tails and of PIPKI. 
PIPKI(S1!1,Y644~U).SRC(SH3!1,Y419-P) -> PIPKKSl ! 1,Y644-P) .SRC(SH3! 1,Y419-P) k8 



Rule 9: Integrin phosphorylation Integrins have to be in an open conformation, SRC bound, but 
not necessarily bound by ligand. 

SRC(SH3! 1,Y419~P) .INT(BD! 1,NPXY~U) -> SRC(SH3! 1,Y419~P) .INT(BD! 1,NPXY~P) k9 



Rule 10: Tal and PIPKI interaction In the quiescent state most talin cannot bind integrins because 
of self-inhibitory interactions (Goksoy et al., 2008). The inhibited forms can, however, still bind PIPKI. 
The affinity between PIPKI and talin has been determined as 170 nM (do Pereda et al., 2005 Barsukov 



et al., 2003). Since PIPKI exists as a dimer, we model two talin binding sites leading to TALiPIPKI 



dimers and trimers. Without loss of generality we assume a sequential binding. 
TAL(Loc-PM,BD) + PIPKI (Loc~PM,Tl,T2) <-> 
TAL(Loc~PM,BD!l) .PIPKI(Loc~PM,Tl!l,T2) klOon, klOoff 
TAL(Loc~PM,BD) + PIPKI (Loc~PM, Tl !+,T2) <-> 
TAL(Loc~PM,BD!l) .PIPKI(Loc~PM,Tl!+,T2!l) klOon, klOoff 
TAL(Loc~CP,BD) + PIPKI (Loc~CP , Tl , T2) <-> 
TAL(Loc~CP,BD!l) .PIPKI(Loc-CP,Tl!l,T2) klOon, klOoff 
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TAL(Loc-CP,BD) + PIPKI(Loc-CP,Tl !+,T2) <-> 
TAL(Loc~CP,BD!l) .PIPKI(Loc-CP,Tl!+,T2!l) klOon, klOoff 



Rule 11a: Activation of talin Most cellular talin is unable to bind integrins because of inhibitory 



self-interactions (Goksoy et al., 2008). Talin can be activated by PIPKI which produces PIP2 (Paolo 



et al., 2002 Ling et al., 2002). PIP2 is highly unstable (Hilgemann(2007), and its production, decay and 



diffusion are not modelled explicitly. We rather assume that talin must be bound to membrane-bound, 
active PIPKI for activation. The rate of PIPKI-dependent PIP2 production and talin activation has not 
been measured, and we use fclla — 0.1s~^. 
TAL(Loc-PM,Conf~C,BD! 1) .PIPKI(Loc~PM,Tl! 1,Y644~P) 
TAL(Loc~PM,Conf~D,BD! 1) .PIPKI(Loc~PM,Tl! 1,Y644-P) 
TAL(Loc~PM,Conf-C,BD! 1) .PIPKI(Loc~PM,T2! 1,Y644-P) 
TAL(Loc~PM,Conf ~D,BD! 1) .PIPKI(Loc~PM,T2! 1,Y644~P) 



-> 

klla 
-> 

klla 



Rule lib: Activation of talin We also assume a PIPKI independent activation of talin at the 
membrane with a very low rate of kllb = 10^'^. 
TAL(Loc~PM,Conf ~C) -> TAL(Loc~PM,Conf -0) kllb 



Rule 12: Inactivation of talin Only if it is unbound 
TAL(Conf~0,BD) -> TAL(Conf ~C,BD) kl2 



Rule 13: Src de-phosphorylation We include constitutive Src, integrin, and PIPKI dephosphoryla- 
tion at rates 0.1 s^^. These rates have not been measured directly in experiments, but data exist that 
provide some bounds. Thus in talin knock-out cells Src activity drops to 50 % of its maximal value within 
15-20 minutes. For large dephosphorylation rates there is little Src activation and the system rapidly 
equilibrates. For small dephosphorylation rates rapid Src deactivation is impossible. The rate of Src 
dephosphorylation must therefore be about fcl3 = 0.1s~^. 
SRC(Conf-0,Y419-P) -> SRC(Conf-0,Y419~U) kl3 



Rule 14: PIPKI de-phosphorylation Must be unbound by Src 
PIPKI(Y644~P,S1) -> PIPKI(Y644~U,S1) kl4 



Rule 15: integrin dephosphorylation Must be unbound by Src 
INT(NPXY~P,BD) -> INT(NPXY~U,BD) kl5 



Rule 16a: talin shuttling 

TAL(Loc~CP,BD) <-> TAL(Loc-PM,BD) kl6on, kl6off 



Rule 16b: PIPKI shuttling incase of dimers 
TAL(Loc-CP,BD! 1) . PIPKI(Loc-CP , Tl ! 1 , T2 , SI) <-> 
TAL(Loc~PM,BD! 1) .PIPKI(Loc~PM,Tl ! 1,T2,S1) kl6on, kl6off 



Rule 16c: DOK shuttling (same parameters as for talin) 
DOK(Loc-CP,BD) <-> DOK(Loc~PM,BD) kl6on, kl6off 
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Rule 17a: PIPKI shuttling Phosphatidylinositol phosphate kinase type I7 (PIPKI) appear not to 
bind well to membrane lipids and PIPKI is therefore assumed to be homogeneously distributed between 
membrane compartment and cytoplasm (Giudici at al., 2006), i.e PIPKI shuttles between membrane and 
cytoplasm at rate 1 s^^ in both directions. 

PIPKI (Loc~CP,Tl,T2, SI) <-> PIPKI (Loc-PM,Tl , T2 , SI) kl6on, kl7offA 



Rule 17b: PIPKI shuttling in case of trimers 

TAL(Loc-CP,BD! 1) . PIPKI (Loc~CP , Tl ! 1,T2!2,S1) .TAL(Loc~CP,BD!2) <-> 
TAL(Loc-PM,BD! 1) . PIPKI (Loc-PM, Tl ! 1 , T2 ! 2 , SI) . TAL(Loc~PM,BD ! 2) kl6on, kl7offB 



Compartment modeling 

We consider two reaction compartments: cytoplasm (CP) and plasma membrane (PM). TAL, PIPKI and 
DOK can shuttle between both compartments. Since the CP volume is about 20 fold larger compared 
with the effective PM reaction volume, the shuttling rates have to account for the change in species 
concentration due to shuttling. Therefore, we normalize the compartment exchange fluxes on a per- 
species basis as detailed below. Note, that this approach differs from the standard SBML formulation, 
which relies on an extensive description of species levels (i.e. species levels are treated as absolute numbers 
and not as concentrations) (Hucka et a/. (2003). 

Let A and B be two compartments with reaction volume Va and Vb, respectively. Let ua be the 
number of particles in compartment A and ns the number of particles in compartment B. For simplicity 
we assume a constant number of particles. 



n-T = nA + riB 



The concentrations in the compartments are denoted as xa = tia/Va and xb = tib/Vb- We want to 
describe the change in concentration due to the flux of particles between the compartments. We assume 
large particle numbers, such that the problem can be formulated in terms of a set of differential equations. 

XA = -OlA^B + a.B^A (1) 

XB = -Pb^a + Pa^b (2) 

We will assume mass-action kinetics, i.e., all fluxes are proportional to the concentrations of the reacting 
species. 



OiA^B = aiXA 

aB^A = a2XB 
Pb^a = biXB 
Pa^B = b2XA 

Due to conservation of particles fiA + tib = 0, and therefore 

Va Vb < 

xaVa + xbVb = -aiUA + 02—^3 - biUB + 62— = 0. 

VB Va 

The last equation can be rearranged to 

n • u = 0, 
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where n ~ {ua^ ns) and v — (62-^ ~ ai, 0.2^ — This equation must hold for ah n. As an example, 

consider the two equations ni ■ v = and n2 ■ v = for linearly independent vectors ni and n2, which 
can only hold true for v = (0, 0). This is therefore the general solution for arbitrary n. It follows 

bi =a2^ 

Inserting these into the original equation for xb leads to 

±B = -{aiVAXB + a2VAXA)/VB. 

Defining the new parameters oi = QiVa and 02 — 02 Va, we can rewrite the original equations in terms 
of the new parameters as 

XA ^ {~diXA + d2XB)/VA (3) 

Xb = {-d2XB + diXA)/VB- (4) 

The dimension of di and 02 are volume over time and are explicitly referencing Va- However, we can by 
the same approach define parameters bi ~ feiVg and 62 = ^2^3 referencing volume Vb and derive the 
equations in terms of bi and 62- 

XA = {-hXA + biXB)/VA 
Xb = i~blXB + b2XA)/VB 

From this it is apparent that oi = 62 = kA-^B and 02 = 61 = fes-fA- 

Particle numbers might not be conserved due to production or decay occurring within compartments. 
However, in mass-action kinetics velocities are proportional to concentrations: changing the reaction 
volume while keeping the concentrations fixed, does not change the velocity of a reaction. Therefore, the 
transformations ai — >■ di and 02 — 02 do not affect other fluxes within the compartments, whose rates 
are volume independent. Therefore the general form of the reaction kinetics including production, decay 
and compartment exchange is: 

XA = {~diXA + d2XB)/VA + fprod{xA) - fdeg{xA) (5) 
XB = {-d2XB + diXA)/VB + fprodixB) - fdeg{xB)- (6) 



Evaluation of sample size 

Exhaustive sampling of the biologically feasible parameter range is computationally infeasible. Since 
most of the parameter combinations lead to qualitatively similar model dynamics we are interested in the 
effective number of samples needed to get a comprehensive picture of the possible model dynamics. To this 
end we monitor the convergence of the time-dependent mean of the model observables for an increasing 
number of samples by a blocking procedure. Samples of size = {10, 10^, 10^, 10"*, 10^} are divided into 
10 blocks. For each sample size we calculate the overall mean as well as the standard deviation of the 
means of all 10 blocks. Next, we compute the coefhcient of variation (CV), i.e., the standard deviation 
divided by the mean, for each sample size. This gives us a CV value for each sample size, time point and 



model observable (Figure SI). The analysis shows, that the average signaling dynamics is sufficiently 
well captured with 10^ samples, i.e. CV < 0.1 for all observables and time points. 



7 



Supplementary References 

Akiyama et a/. (1985. Akiyama SK, Hasegawa E, Hasegawa T, Yamada KM (1985) The interaction of 
fibronectin fragments with fibroblastic ceUs. J Biol Chem 260: 13256-60 

Anthis et aL(2009. Anthis NJ, Haling JR, Oxlcy CL, Memo M, Wegener KL, Lim CJ, Ginsberg MH, 
Campbell ID (2009) Beta integrin tyrosine phosphorylation is a conserved mechanism for regulating 
talin-induced integrin activation. J Biol Chem 284: 36700-36710 

Arias-Salgado et a/.(2003. Arias-Salgado EG, Lizano S, Sarkar S, Brugge JS, Ginsberg MH, Shattil SJ 
(2003) Src kinase activation by direct interaction with the integrin beta cytoplasmic domain. Proc 
Natl Acad Sci USA 100: 13298-13302 

Barsukov et aL(2003. Barsukov IL, Prescot A, Bate N, Patel B, Floyd DN, Bhanji N, Bagshaw CR, 
Letinic K, Di Paolo G, De Camilli P, Roberts GCK, Critchley DR (2003) Phosphatidylinositol 
phosphate kinase type Igamma and betal-intcgrin cytoplasmic domain bind to the same region in 
the talin PERM domain. J Biol Chem 278: 31202-9 

Ben(1978. Bell G (1978) Models for the specific adhesion of cells to cells. Science 200: 618-627 

Calderwood et aZ.(2002. Calderwood DA, Yan B, de Pereda JM, Alvarez BG, Fujioka Y, Liddington 
RC, Ginsberg MH (2002) The phosphotyrosine binding-like domain of talin activates integrins. J 
Biol Chem 277: 21749-21758 

Cox et a/. (2001. Cox EA, Sastry SK, Huttenlocher A (2001) Integrin-mediated adhesion regulates cell 
polarity and membrane protrusion through the Rho family of GTPases. Molecular Biology of the 
Cell 12: 265-77 

de Pereda et aL(2005. de Pereda JM, Wegener KL, Santelh E, Bate N, Ginsberg MH, Critchley DR, 
Campbell ID, Liddington RC (2005) Structural basis for phosphatidylinositol phosphate kinase 
type Igamma binding to talin at focal adhesions. J Biol Chem 280: 8381-8386 

Faull & Ginsbcrg(1995. FauU RJ, Ginsberg MH (1995) Dynamic regulation of integrins. Stem Cells 
13: 38 46 

Giudici et aL(2006. Giudici ML, Lee K, Lim R, Irvine RF (2006) The intracellular localisation and 
mobility of Type Igamma phosphatidylinositol 4P 5-kinase splice variants. FEBS Lett 580: 6933- 
7 

Goksoy et a/. (2008. Goksoy E, Ma YQ, Wang X, Kong X, Perera D, Plow EE, Qin J (2008) Structural 
basis for the autoinhibition of talin in regulating integrin activation. Mol Cell 31: 124-33 

Harrison(2003. Harrison SC (2003) Variation on an Src-like theme. Cell 112: 737-40 

Hilgemann(2007. Hilgemann DW (2007) Local PIP (2) signals: when, where, and how? Pflugers Arch 
455: 55-67 

Hucka et aL(2003. Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein 
BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin 
II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, et al. (2003) The systems biology markup 
language (SBML): a medium for representation and exchange of biochemical network models. 
Bioinformatics 19: 524-31 

Hynes et al.{2002. Hynes RO, Lively JC, McCarty JH, Taverna D, Francis SE, Hodivala-Dilke K, Xiao 
Q (2002) The diverse roles of integrins and their ligands in angiogenesis. Cold Spring Harb Symp 
Quant Biol 67: 143-53 



8 



Ibcr & Campbcll(2006. Iber D, Campbell ID (2006) Integrin activation-the importance of a positive 
feedback. Bull Math Biol 68: 945 56 

Kuo & Lauffenburgcr(1993. Kuo SC, Lauffenburgcr DA (1993) Relationship between receptor/ligand 
binding affinity and adhesion strength. Biophys J 65: 2191 200 

Ling et a/. (2002. Ling K, Doughman RL, Firestone A J, Bunce MW, Anderson RA (2002) Type I 
gamma phosphatidylinositol phosphate kinase targets and regulates focal adhesions. Nature 420: 
89-93 

Ling et aL(2003. Ling K, Doughman RL, Iyer W, Firestone AJ, Bairstow SF, Mosher DF, Schaller 
MD, Anderson RA (2003) Tyrosine phosphorylation of type Igamma phosphatidylinositol phos- 
pliate kinase by Src regulates an integrin-talin switch. J Cell Biol 163: 1339-1349 

Moser et a/. (2009. Moser M, Legate KR, Zent R, Fassler R (2009) The tail of integrins, talin, and 
kindlins. Science 324: 895-899 

Moy et a/. (1999. Moy VT, Jiao Y, Hillmann T, Lehmann H, Sano T (1999) Adhesion energy of 
receptor-mediated interaction measured by elastic deformation. Biophys J 76: 1632-1638 

Neff et a/.(1982. Neff NT, Lowrcy C, Decker C, Tovar A, Danisky C, Buck C. Horwitz AF (1982) A 
monoclonal antibody detaches embryonic skeletal muscle from extracellular matrices. J Cell Biol 
95: 654-66 

Oxley et a/. (2008. Oxley CL, Anthis NJ, Lowe ED, Vakonakis I, Campbell ID, Wegener KL (2008) 
An integrin phosphorylation switch: the effect of beta3 integrin tail phosphorylation on Dokl and 
talin binding. J Biol Chem 283: 5420-5426 

Paolo et a/. (2002. Paolo CD, Pellegrini L, Lctinic K, Cestra G, Zoncu R, Voronov S, Chang S, Guo 
J, Wenk MR, Camilli PD (2002) Recruitment and regulation of phosphatidylinositol phosphate 
kinase type 1 gamma by the FERM domain of talin. Nature 420: 85-9 

Plow & Ginsberg(1981. Plow EF, Ginsberg MH (1981) Specific and saturable binding of plasma fi- 
bronectin to thrombin-stimulated human platelets. J Biol Chem 256: 9477-82 

Roskoski(2004. Roskoski R Jr (2004) Src protein-tyrosine kinase structure and regulation. Biochem 
Biophys Res Commun 324: 1155-64 

Saltel et a?. (2009. Saltel F, Morticr E, Hytoncn VP, Jacquier MC, Zimmermann P, Vogel V, Liu W, 
Wehrle-Haller B (2009) New PI(4,5)P2- and membrane proximal integrin-binding motifs in the 
talin head control beta3-integrin clustering. J Cell Biol 187: 715-731 

Songyang et aZ.(1995. Songyang Z, Carraway KL 3rd, Eck MJ, Harrison SC, Feldman RA, Mohammadi 

M, Schlcssinger J, Hubbard SR, Smith DP, Eng C (1995) Cataljd:ic specificity of protein-tyrosine 
kinases is critical for selective signalling. Nature 373: 536-9 

Tadokoro et aL(2003. Tadokoro S, Shattil SJ, Eto K, Tai V, Liddington RC, de Pereda JM, Ginsberg 
MH, Calderwood DA (2003) Talin binding to integrin beta tails: a final common step in integrin 
activation. Science 302: 103-6 

Vitte et al.(2004. Vitte J, Benoliel AM, Eymeric P, Bongrand P, Pierres A (2004) Beta-1 integrin- 
mediated adhesion may be initiated by multiple incomplete bonds, thus accounting for the func- 
tional importance of receptor clustering. Biophys J 86: 4059-74 



9 



Wiseman et oL(2004. Wiseman PW, Brown CM, Webb DJ, Hcbcrt B, Johnson NL, Squier JA, Ellis- 
man MH, Horwitz AF (2004) Spatial mapping of integrin interactions and dynamics during cell 
migration by image correlation microscopy. J Cell Sci 117: 5521-34 



10 



Open Integrins 



Phosphorylated Integrins 




1 sec 1 min 
Time 



100000 




1 sec 1 min 1 h 
Time 



Ligand-bound Integrins 



Dok-bound Integrins 




1 sec 1 min 
Time 



2.5 










100000 


2 








CD 


10000 




N 




1.5 


in 






CD 
CL 


1000 


1 


E 




n3 






cn 


100 


0.5 










10 




1 sec 1 min 
Time 



Talin-bound Integrins 



Phosphorylated SRC 




1 sec 1 min 
Time 



100000 




1 sec 1 min 
Time 



Figure SI. Mean signaling dynamics in dependence of the sample size. Each plot shows the 
coefficient of variation (standard deviation/mean) in dcpendance of time (horizontal axis) and sample 
size ( vertical axis) for a selected observable. Mean and standard deviation are calculated by a blocking 
procedure as described in the text. For a sample size of 10^ the coefficient of variation is below 0.1 % 
for all time points and observations. 
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Table SI. Mean parameter values of group 1 (high TAL:INT, Figure 3 main text) and group 2 (high 
DOKiINT, Figure 3 main text). On-rates are related to the KD value mentioned in the main text as 
kon = koff/KD. 



